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Abstract 


Accurate flux density calibration is essential for precise analysis and interpretation of observations across different 
observation modes and instruments. In this research, we first introduce the flux calibration model that incorporated 
in HiFAST pipeline, and designed for processing HI 21cm spectra. Furthermore, we investigate different 
calibration techniques and assess the dependence of the gain parameter on the time and environmental factors. A 
comparison is carried out in various observation modes (e.g., tracking and scanning modes) to determine the flux 
density gain (G), revealing insignificant discrepancies in G among different methods. Long-term monitoring data 
shows a linear correlation between G and atmospheric temperature. After subtracting the G—Temperature 
dependence, the dispersion of G is reduced to <3% over a one-year timescale. The stability of the receiver 
response of Five-hundred-meter Aperture Spherical radio Telescope (FAST) is considered sufficient to facilitate 
HI observations that can accommodate a moderate error in flux calibration (e.g., > ~5%) when utilizing a 
constant G for calibration purposes. Our study will serve as a useful addition to the results provided by Jiang et al. 
Detailed measurement of G for the 19 beams of FAST, covering the frequency range 1000-1500 MHz, can be 
found on the HiFAST homepage: https: / /hifast.readthedocs.io /fluxgain. 
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1. Introduction 


Observational data acquired from radio telescopes are often 
subject to a multitude of origins of biases, including weather 
conditions, radio frequency interference (RFI), hardware drift, 
deformation of the telescope panels, etc. Hence, calibrating the 
acquired data is imperative to derive reliable results. Well- 
established methods utilize robust and stable radio sources as 
calibrators, with extensive literature documenting their efficacy 
(Kuzmin 2013; Baars 1973; Baars et al. 1977; Boothroyd et al. 
2011; Salter 2000). Fluctuations in the hardware configuration 
of a system can induce variations in the gain of a radio 
telescope, necessitating regular calibration of the electronic 
system (Maddalena 2004). The choice of a calibration method 
is contingent upon several factors, including the observing 
strategy (e.g., observation modes, data sampling time), 
frequency range, and telescope performance. Hence, selecting 
the appropriate calibration method is pivotal for obtaining 
accurate results in radio astronomy observations. 


l Corresponding authors. 


Among the various calibration processes, precise flux density 
calibration is particularly crucial, especially for highly sensitive 
single-dish telescopes. Single-dish telescopes serve as indis- 
pensable tools for probing the large-scale, low-density distribu- 
tion of gas (van Zee et al. 1997). The Five-hundred-meter 
Aperture Spherical radio Telescope (FAST) completed its 
commissioning phase in 2019 (Nan 2006; Nan et al. 2011; Jiang 
et al. 2019). Featuring a 19-beam receiver design and 
unparalleled sensitivity, FAST stands as an optimal instrument 
for conducting large-scale HI surveys (Li et al. 2018; Duffy 
et al. 2008). Several forecasts regarding FAST HI spectral 
research have been articulated, including its capability to 
spearhead the next generation of HI blind sky surveys aimed 
at detecting high-redshift galaxies and other faint radio sources 
(Duffy et al. 2008). Ongoing initiatives such as the Commensal 
Radio Astronomy FAST Survey (CRAFTS, Li et al. 2018) are 
underway, with Zhang et al. (2019) providing insights into the 
survey’s extragalactic HI sensitivity and projecting FAST’s 
future survey capabilities. By combining observational data from 
FAST with that from WSRT and VLA, Wang et al. (2023, 2024) 
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have achieved a more comprehensive understanding of the 
diffuse H I distribution surrounding interacting galaxies. 
Additionally, the FASHI source catalog published by Zhang 
et al. (2023) highlights FAST's exceptional performance in blind 
survey projects. 

To realize the scientific objectives outlined above, precise flux 
density calibration is imperative for the sky surveys and 
observations conducted using FAST. However, the existing data 
open-source data-processing pipelines for FAST lack a systematic 
flux density calibration procedure. Therefore, we have developed 
and integrated a module into the HiFAST pipeline specifically 
designed for flux density calibration of observational data, 
followed by an evaluation of its performance. 

Several observation modes are offered to observe calibrators 
by FAST, including OnOff, MultiBeamCalibration, and Drift. 
Selecting the most suitable calibration mode for a given 
observation is crucial. Furthermore, the stability and variability 
of the gain depend on both the frequency and the directionality 
of the telescope. As FAST utilizes a transformable paraboloid 
surface for observation, the mirror undergoes deformation at 
varying elevations and azimuth angles, potentially leading to 
fluctuations in observation gain throughout the observation 
process. The stability of the FAST electronic system meets the 
accuracy requirement (<1%) for observations with calibrators 
(Jiang et al. 2020). Despite the importance of calibration, its 
implementation is often constrained by limited observation 
time and the availability of calibrators. Notably, calibrating 
each of the 19 beams individually consumes a significant 
amount of observation time. Thus, in scenarios where 
calibration is unfeasible, a fixed G value is adopted for flux 
density calibration. Jiang et al. (2020) provided a version of the 
FAST G table (hereafter G 2020) with a fitting function of zenith 
angle (Oza) at specific frequencies. It is anticipated that 
observation modes and external environmental factors might 
also affect G. To date, the accuracy of flux calibration across 
various observation modes has not been thoroughly investi- 
gated. Furthermore, there has been no examination of the long- 
term stability of the gain parameter of FAST. This study aims 
to investigate different calibration techniques and assess the 
accuracy of the gain parameter over an extended duration. This 
might serve as an addition to G 020- 

In this paper, the calibration selection criteria and observation 
strategy are delineated in Section 2. Following that, the data 
reduction procedures for various observation modes used in 
calibration are expounded upon in Section 3. Section 4 presents 
a comparative analysis of the G values obtained from different 
observation modes and examines the gain curve across the 19 
beams. Additionally, the section delves into the longitudinal 
variations of G observed during extended observations. After 
this, Section 5 encompasses in-depth discussions on the 
outcomes. Finally, Section 6 offers a succinct summary. 


Liu et al. 


2. Calibration Data 


In this section, we will first introduce the choice of calibrator 
sources, then the observation modes often used for calibration. 
Finally, the introduction to the calibration data used in this 
paper will be presented. 


2.1. The Choice of Flux Calibrators 


The choice of absolute flux density calibrators is based on four 
criteria: (i) the source should have been regularly monitored to 
ensure that their L-band fluxes do not vary with time, or within a 
limited variation range (Perley & Butler 2017); (ii) spatial 
confusion in calibration should be avoided, meaning that there 
should be no other continuum sources listed in the NVSS catalog 
provided by Condon et al. (1998) close to the calibrator; (iii) the 
coordinates of the calibrator should be as close as possible to the 
target sources; and (iv) the target should be a point source, with its 
angular extent much less («K 10% or less) than the beam size. 

For our test sky coverage around the M31 galaxy, 3C48 (also 
known as J0O137+3309) and 3C286 (also known as J1331 
+3030), were chosen as calibrators. These sources have been 
extensively observed by VLA, VLBA, MERLIN, and EVN 
(An et al. 2010) and are usually used as calibrators in 
observations with VLA and other telescopes due to their 
stability flux and small angular size (Perley & Butler 2013a). 
The flux densities of 3C48 and 3C286 can be accurately fitted 
by a polynomial function: 


log[S.(v)] = ao + alog) + arflog()P + as[log()P , 
(1) 


where S denotes the flux density in Jy at the frequency v in 
GHz. The coefficients used are dg = 1.3253, a; = —0.7553, 
a = —0.1914, and a3 = 0.0498 for 3C48, and ap = 1.2481, 
a, = —0.4507, az = —0.1798, and a3 = 0.0357 for 3C286 from 
Perley & Butler (2017). In the frequency range of 1000-1500 MHz, 
the flux density of 3C48 exhibits negligible polarization, with 
a polarization difference of 0.3% at 1050 MHz and 0.5% 
at 1450 MHz as reported by Perley & Butler (2013b). In 
contrast, 3C286 serves as a reliable polarization calibrator due 
to its stable polarization degrees, angles, and total flux density 
within this frequency range (Perley & Butler 2013a). 

The two calibrator sources that we selected have very similar 
values of zenith angles in culmination time: ~8° for 3C48, and 
~5° for 3C286. The difference is less than 3°. The observations 
of the calibrators are always processed near their culmination 
time, and the changes on zenith angles during observation are 
always less than 1°2 for any DriftWithAngle, OnOff, Multi- 
beamCalibration, or MultibeamOTF mode observations. Accord- 
ing to G pozo, this difference has minimal impact on the gain. So 
no additional correction on the G during these observations is 
needed. Anyway, for all mode observation data, we followed the 
correction table provided in Gy929 and corrected the gain 
parameter in all of our samples to the zenith angle of 10.15°. 
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Figure 1. The schematic of observation modes used in calibration: 
DriftWithAngle (left-top), MultiBeamOTF (left-bottom), MultiBeamCalibra- 
tion (right-top), OnOff (right-bottom). The geometry of 19 beams is shown in 
blue dot (start position) and orange dot (end position), the red cross is the 
position of the calibrator, and the blue dashed lines show the desired traces of 
Beam 1. The relative positions of M01 and M08 are also marked on the first 
panel. 


2.2. The Modes of Calibrations 


The accuracy of the calibration highly relies on the calibration 
mode, i.e., the observing strategy of the calibrators. The 
calibration process first involves a pointing calibration. If the 
calibrator source is significantly smaller than the beam size, such 
that it can be treated as a point source, the flux variation can be 
effectively represented by the beam response of the telescope 
(O’Neil 2002). The beam response of the FAST receivers 
exhibits high circular symmetry and can be described using a 
two-dimensional Gaussian model. Then the flux density of the 
calibrators can be determined by fitting a Gaussian function to 
the observed response changes in a one-dimensional scan, as 
demonstrated in Jiang et al. (2020) and Xi et al. (2022). 

FAST provides various observation modes for calibration, 
and we examined four of them: DriftWithAngle, Multi- 
BeamOTF, MultiBeamCalibration, and OnOff. Figure 1 
provides a brief illustration of these four observation modes. 
The desired trajectory of Beam 1 is denoted by blue dashed 
lines, while the configuration of the 19 beams is represented by 
blue and orange dots. In the following, we describe the four 
observing modes of calibrations adopted by FAST: 

OnOff. This mode enables a rapid positional switch between 
the source and reference positions with a slewing time of 30 s. 
Both positions are observed with the same beam and with the 
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same duration. For multiple beam receivers, such as the L-band 
19-beam receiver on board FAST, the separation between 
positions can be optimized for the simultaneous measurement of 
two beams. In our test, beams | (the center beam) and 8 (one of 
the outermost beams) were chosen for calibration, and the 
duration time was set to 60 s to reduce the error caused by the 
telescope’s morphology instability at the start of the observation. 
A high noise diode signal with a synchronized period of 
2.01326592 s was injected during the observation, with a signal 
strength of approximately 11 K. In total, it took 150 s to calibrate 
the two beams, including 60 s of tracking on the source, 60 s of 
tracking off the source, and 30 s of slewing time. 

MultiBeamCalibration. This mode is able to calibrate all 19 
beams automatically by tracking the source with each beam in 
succession, following the order: 1—1—2—3-4—5-6-7-19-8 
9-10-11-12—13-—14—15-16-17-18. Each beam requires 40 s 
to switch, and the same amount of time is needed to track the 
beam and set the noise diodes. This process takes a total of 
1960 s to complete. 

DriftWithAngle. In this mode, the receiver was rotated at an 
angle of 0° and the calibrator was scanned at a rate of 15” s7! 
with scan lines of length 67/5. This setup enabled the centers of 
beams 1, 2, 5, 8, and 14 to pass the calibrator in one scan with a 
sufficient baseline on either side. The entire observation took 
270 s to calibrate the five beams and the high noise diode signal 
was injected only at the start and end of the scanning lines with 
a duration of 1.00663296 s. 

MultiBeamOTF. Similar to the DriftWithAngle mode, this 
mode repeats five times with a rotation angle of 0° to calibrate all 
19 beams. Scanning lines are spaced at an interval of 4/97, and 
the same noise diode injection technique as the DriftWithAngle 
method is used. The total time needed for one MultiBeamOTF 
observation to calibrate all 19 beams is 1260 s. 

The calibrator observations for all modes were recorded using 
the spectral backend with a sampling rate of 0.201326592 s and 
a channel width of 7.62939 kHz. 


2.3. Pointing Accuracy 


We should be aware that the telescope’s pointing accuracy can 
also influence the measurement of the parameter G. The pointing 
error of FAST included errors from both the surface precision of 
the reflector system and the accuracy of the feed support system. 
The measurement error of the reflector system is primarily 
influenced by spherical aberration, deviation in measurement 
shafting, atmospheric refraction, and the accuracy of the 
measurement reference net and the measurement errors of the 
feed support system come mainly from the error of the 
measurement reference net, the target error of the total station, 
the random error of the total station, and the time delay (Jiang 
et al. 2019). Raster scan observations have been conducted to 
measure the pointing errors of FAST through raster scan 
observations (Jiang et al. 2020). Jiang et al. (2020) used a 2-D 
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Figure 2. Relative position offset measured in OnOff mode observations. Blue 
plus points show the coordinates of sampling points in the source area observed 
on 2021 September 9, distinguished by the observation time record by FAST. 
The red shadow area shows the region that is included in 16” of 3C48. The red 
star presents the center of the calibrator, 3C48. 


Gaussian model to model the pointing errors. This model is 
suitable because the central beam does not exhibit significant 
beam ellipticity or coma features. The surface precision of the 
FAST reflector system is expected to reach an rms of 8” (Jiang 
et al. 2019). Jiang et al. (2020) confirmed that the pointing error 
of the FAST beams is smaller than 16”, with an rms of 7”9. An 
example of the coordinate distribution of the sample points of the 
OnOff observation on 2021 September 9, is shown in Figure 2. 
The red shadowed area indicates the filtered range of 16”. The 
pointing remains within the telescope’s pointing error range 
(16”). However, occasionally, a few points may exceed this 
range. To reduce those effects, we filtered out these outlier 
sample points during data processing. Through the analysis of 
scanning calibration conducted at various frequencies, the 
fluctuation of gain values near the beam’s center with pointing 
offset has been investigated. This examination reveals the impact 
of pointing accuracy on calibration observations in 1-dimension, 
as depicted in Figure 3. We found that the flux attenuation caused 
by the FAST pointing error was always less than 1%, so this error 
was not taken into account in our subsequent analysis. 

The flux bias of pointing rms errors, which occur when 
utilizing a single-dish telescope to observe a point-like source, 
is corrected by a correction coefficient fg derived from the 
following equation (Papadopoulos et al. 2010): 


F. 2 
=) (2) 


,V)=1+4 8ln2 
folo ni j [st 
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Figure 3. Gaussian function fitting of the 1-dimensional scan of MO1 at a 
different frequency. The x—axis is the R.A. offset of sample points with the 
calibrator or beam center, and the y-axis shows the difference of detected 
power in percentage. 


where ø is the pointing rms error of FAST which is 8”, we use 
the beam size(Oyppw(v)) of FAST in Table 2 of Jiang et al. 
(2020). The value of fg is 1.0079 at 1100 MHz and 1.0095 at 
1400 MHz. 


2.4. Observation Data 


The calibrator’s observational data serves two research 
purposes: evaluating various calibration methods and examin- 
ing gain stability over extended periods. For the first part, we 
performed continuous observations of 3C48 using three 
different observation modes to assess the flux density 
calibration accuracy under varying modes. These observations 
were conducted on 2021 September 9, and 2022 June 29, 
respectively. During the 2021 September 9, observations, we 
utilized the MultiBeamCalibration, DriftWithAngle, and OnOff 
modes, while during the 2022 June 29, observations, we 
replaced the DriftWithAngle mode with the MultiBeamOTF 
mode to calibrate all beams of FAST. For the second part, the 
long-term calibrator observations for the M31 halo HI survey 
were carried out subsequent to each DriftWithAngle survey 
observation. We performed flux calibration for each beam 
roughly every 10 days, accumulating the calibration data over a 
period of more than 175 days (for the center beam, also M01), 
from 2020 September to 2022 September. Additionally, a 
limited amount of publicly accessible 3C286 data was 
employed to enhance our study when 3C48 was not proper 
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for observations, with these observations being carried out in 
MultibeamCalibration mode. 


3. Data Reduction 


All data used in this paper are processed with our pipeline 
HiFAST (Jing et al. 2024), including frequency-dependent 
noise diode calibration, baseline fitting, standing wave 
removal, flux density calibration, stray radiation correction, 
and gridding to produce data cubes. For further details, please 
consult the series paper. This section is dedicated solely to the 
model for flux density calibration. 

The HiFAST.FluxGain module of the HiFAST pipeline 
is used to calibrate the flux density gain of FAST. First, the 
observed power of the calibrator (Pobs) is converted to the 
antenna temperature (T,,) by using the noise diode signal with a 
fixed known temperature of fae 


Ros) 
cal v) 


T) = 


x Tai), (3) 


where Peai is the recorded power of the injected noise diode. 
With T, obtained for each spectrum per unit frequency, for a 
source with given flux S.(v), G can be expressed as 


PHT w) 
= 4 
Gv) Gh (4) 


TSPC 


where is the antenna temperature of the spectrum in the 


source area, T"F is the spectrum in the reference area. We 
briefly discuss the techniques to obtain Te and 7 
separately in the following two subsections. 


3.1. Tracking Modes 


When using tracking modes such as OnOff and Multi- 
BeamCalibration, the FAST backend records the sampling time 
of the spectrum, which is divided into source (SRC) and 
reference (REF) groups. The mean spectrum in each group is 
used to calculate the antenna temperature of the source and 
reference (TSRS and TR"). By introducing these two values 
into Equation (4), the gain can be determined. 


3.2. Scanning Modes 


In scanning modes such as DriftWithAngle and Multi- 
BeamOTF, the power of the observed spectra will vary 
significantly as the beam passes over the source. This is the 
result of the convolution of the calibrator source and the shape 
of the beam. Generally, the main lobe of the FAST beams can 
be approximated as a Gaussian function (Jiang et al. 2020; Xi 
et al. 2022). To obtain TSPC and TR®", we use the following 


8 https: / /fast.bao.ac.cn/cms/category /telescope_performence_en/noise_ 
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Figure 4. An example of Gaussian function fitting results of DriftWithAngle 
observation. The blue dots represent the sample points, and the red line is the 
best-fit curve. The x—axis is the arc second offset between the sample point and 
the calibrator, and the y—axis is T,. The top panel is fitted with a Gaussian 
function, and the bottom panel is fitted with a linear function. 


function to fit the observed data: 


_ (d-d) 


) 
f(d, v) = TE®(d, v) x e zod + TR (d, v) (5) 


where d is the angular distance offset, dọ is the position of 
calibrator position, o is the standard deviation depended on 
frequency. 

The upper panel of Figure 4 displays an example of a 
Gaussian fit to the temperature variance of Beam 1. When the 
total residual of the fitting reaches its minimum, rs and TREF 
can be determined simultaneously. Assuming the main lobe 
shape follows a Gaussian function may introduce an error of 
approximately 3%, given its significantly more complex nature. 
To reduce this error, we fit separately for TS®© and TR®®, For 
TEE , a linear function is used to fit the reference area, as 
illustrated in the lower part of Figure 4. Then, the equation is 
applied to the data points passing through the calibration source 
to determine TSPC, With TR"* and TR®", we then obtain G with 
Equation (4). 


4. Results 


4.1. Distinction of Gain with Different Observation 
Modes 


To investigate the influence of observation modes on 
calibration accuracy, we conducted continuous observations 
of 3C48 using three modes on two different dates: 2021 
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Figure 5. Gain parameter G of beam 1 among different methods. Symbols are 
for different modes respectively. The results from two different days (2021 
September 9, 2021, and 2022 July 29) are presented by solid blue lines and 
dashed green lines respectively. The red dashed line and error bar are the 
results from Gyj2929 which is observed with the OnOff mode. The gray- 
shadowed region shows the high probability frequency band (>50%) in Zhang 
et al. (2022). 


September 9, (MultiBeamCalibration, DriftWithAngle, and 
OnOff) and 2022 June 29, (replacing DriftWithAngle with 
MultiBeamOTF to obtain G for all beams). The electronic gain 
fluctuations of the FAST system are known to be better than 
1% over hours of observation (Jiang et al. 2020). Additionally, 
the aperture efficiency of FAST varies as a function of 0z4, we 
corrected for the influence of 07, in processing the flux density 
gain data using the fitting results of Jiang et al. (2020). 
Therefore, the observation mode will not significantly affect the 
value of G if the difference in G obtained on the same date is in 
a similar order. 

Figure 5 illustrates the variation in G for Beam 1 across three 
observation modes and two distinct days. The left and right 
y-axes represent the gain parameter and the difference in ratio, 
respectively. Additionally, for comparison purposes, the results 
of Gjoo29 (Jiang et al. 2020) are depicted as red dotted-dashed 
lines. Different line styles (solid and dashed) denote observa- 
tions conducted on different dates, while distinct line colors 
differentiate between various observation modes. 

First, the disparity between the various observation modes is 
quite minimal, with a difference of less than two percent across 
the entire frequency range. Second, the MultiBeamOTF 
observation on 2021 September 9 at 1200 MHz yielded a 
higher G than the other two modes, which is likely due to the 
presence of strong RFI signals mixed in with the observational 
signal. In the shadow range, strong RFI is likely caused by 
communication and navigation satellites (Wang et al. 2021; 
Zhang et al. 2022). So the measurement of the response of 
FAST received at 1000-1500 MHz does not depend on the 
measured method, but could be contaminated by the strong 
RFI. Third, G has a large variation with the observing date. For 
instance, the measurement of G on 2022 June 29 exhibited an 
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increase of nearly 5% compared to the measurement on 2021 
September 9, and a 7% increase compared to the value reported 
in Gj2929. This indicates that the G parameter seems to vary on 
longer timescales of days or even longer. 


4.2. Variation of Gain in Long-term Observation 


To better understand the long-term fluctuations of G, we 
conducted tests on Beam | at around 1400 MHz (averaged with 
a frequency width of 1 MHz) and plotted the results in the left 
panel of Figure 6. The x-axis represents the observation time, 
while the y—axis represents the value of G. We used 3C48 as 
the main calibrator because its coordinates are close to our 
survey area. However, most calibrations in May and June were 
performed with 3C286, as 3C48 was in the Solar avoid- 
ance zone. 

We find that the maximum fluctuation of G is around 8% in 
2022 September compared to the mean value of G. Our 
observations show that the values of G are generally lower in 
winter and higher in late summer, as shown in the right panel of 
Figure 6, we arranged the data by month. This suggests that, in 
addition to system hardware configuration changes, the FAST 
telescope’s response may be influenced by yearly periodical 
factors, such as temperature and air pressure. Next, we roughly 
obtained the relationship between G and observation time using 
polynomial regression, as shown in Figure 6. The polynomial is 
presented below: 


G(t) = ast? + agt* + a3t? + ant? + at + ao (6) 
where t= MJDop; — 59115 is MJD time of observation, the 
parameters of the polynomial are as —=—1.242 x 107", 


a4 = 2.298 x 10°°, a3 =—1.481 x 1076, a= 3.897 x 1074, 
ait = —3.710 x 107°, ay = 16.59. 

To account for the seasonality of the parameter G, we 
examined the atmospheric temperature, pressure, and humidity 
during the calibration source observations. The results indicate 
that the variable G is not sensitive to changes in pressure and 
humidity, while G exhibits a significant relationship with 
atmospheric temperature. Figure 7 shows G as a function of 
atmospheric temperature T, at the observing frequency of 1400 
MHz. We performed a least square linear regression between G 
and atmospheric temperature, with its slope linearly correlated 
with observing frequency, as depicted in Figure 8. Both 
relationships are represented by: 


Goor(v) = G(v) + sv) AT (7) 


where s(v) = 4.7 x 10v —1.4 x 10°, v represents the fre- 
quency of G, and AT denotes the temperature difference 
between the observed value and the standard value. 

After correcting the dependence of the gain G on atmo- 
spheric temperature, we performed a statistical analysis of the 
variance of G over a full calendar year. By setting an 
observation interval window ranging from 1 to 365 days, and 
using a rolling window to traverse all data results, we can 
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Figure 6. Left: Gain parameter G of Beam 1 at 1400 MHz from 2020 September to 2022 September. The diamond marker is G of 3C48 and the square marker is 
3C286. The black dashed line represents the temporal variation of G fitted by polynomial regression. Right: G in the left panel rearranged by months, y-axis is the 


difference between G with mean value in percentage. 
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Figure 7. The correlation between G and atmospheric temperature. The x-axis 
represents temperature and the y-axis represents G. The data points represent G 
of the calibrators at 1400 MHz, and the blue line represents the correlation 
curve fitted using the least squares method. 


calculate the average variance of G over this time interval. This 
method enables us to determine the magnitude of the variation 
in gain values for different observation time intervals. Figure 9 
displays the dispersion of G for different time intervals. The 
solid line represents the 30 dispersion of the original G values 
for each interval, while the dashed line represents G after 
correction using Equation (7). Overall, the dispersion of G 
increases rapidly with longer time intervals. The dispersion of 
the original G exceeds 10% when the time interval is longer 
than 300 days, whereas the corrected G eventually stabilizes 
and remains below 5%, even for intervals close to one year. In 


1050 1100 1150 1200 1250 1300 1350 1400 1450 
Frequency [MHz] 


Figure 8. The relationship between the frequency and slope of the G-T 
relationship. The x-axis represents frequency values, and the y-axis represents 
the values of the slope. Although there are some fluctuations in the variation of 
the slope, overall it still exhibits a linear change with frequency. 


summary, the gain G exhibits remarkable stability, considering 
the systematic errors from the telescope and receiver, as well as 
those introduced during the flux calibration process, which are 
accounted for in the scatter statistics. 

To examine the variability of the scatter for shorter time 
intervals, a zoomed-in version of the dispersion values of G for 
intervals less than 30 days is shown in the inset plot at the 
bottom right of Figure 9. The correction of the G-T relationship 
does not have a significant impact within these shorter time 
intervals. The dispersion of G is approximately 3% over a 
10 days time interval in both situations. This suggests that, in 
the absence of any calibrator observations, the flux calibration 
process could introduce an error of up to 10%. However, this 
error can be reduced to 3% if at least one calibrator observation 
is carried out within 10 days. Moreover, if calibrator 
observations are performed on a daily routine, the error can 
be further reduced to within 2%. On the other hand, if there is a 
long time lag between the calibrator data and the target field 
data, it becomes necessary to correct the gain using the G-T 
relationship. 
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Figure 9. Dispersion of G for different calibrator observation intervals at 1400 
MHz. The dispersion of G represents the average variance of G within a fixed 
interval when traversing all calibrator data using a rolling window approach 
with different observation intervals. The solid line represents the 3a dispersion 
of the original G value for each interval in percentage, while the dashed line 
indicates the G corrected by the G-T relationship. The inset plot at the bottom 
right shows the dispersion of G within the calibration interval of 30 days. 


4.3. Gain Curve Variation Among the Multiple Beams 


During our HI deep survey around the M31 halo region, we 
performed flux calibration for each beam roughly every 10 
days, accumulating the calibration data over a period of more 
than 175 days (for the center beam, also MO1), from 2020 
September to 2022 September. Leveraging these observations, 
we were able to determine the frequency dependency of the G 
(referred to as the gain curve) for every beam. Figure 10 shows 
the gain curve for M01. The solid line represents the mean 
value of G at each frequency channel, while the blue shadowed 
region signifies the dispersion of G post-adjustment. The 
results for Gyo929 are represented by red symbols with error 
bars. In general, the gain parameter gradually declines as the 
frequency surges, except for a small hump around a frequency 
of 1150 MHz. Our findings align with Gj2929 within the 1 
sigma error of the data adjusted by G-T relation. Except for a 
single data point, again at a frequency of 1150 MHz. This 
aberration could be ascribed to considerable RFI contamina- 
tion. The likelihood of RFI emerging around this frequency 
point is roughly 80% and 50% for the XX and YY channels, 
respectively, as indicated by (Zhang et al. 2022). 

Furthermore, in Figure 11, we elucidate the probability 
density function (PDF) of the fractional difference of the G 
parameter around 1400 MHz normalized by the average value, 
in percentage terms. The results with and without G—T 
adjustments are denoted by solid and dashed histograms, 
respectively. For the corrected results, one sigma is merely 
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Figure 10. Frequency dependence of G of Beam 1. The blue line is the mean 
value of flux density gain in all calibration data corrected from the G-T 
relationship, and the pale region shows the 1ø range of corrected—G fluctuation. 
The red dot and error bar are G2929. The right y-axis shows the difference of G 
in percentage. The gray-shadowed region shows the high probability frequency 
band (>50%) in Zhang et al. (2022). 
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Figure 11. The probability distribution function (PDF) of parameter G 
normalized by the mean value around 1400 Mhz. The dashed blue line 
collected original data from every calibration, while the solid line shows data 
corrected by the G-T relationship. The 1ø ranges for the two cases are 3.2% 
and 2.15% respectively. 


about 2.15%, denoted by the shadowed area. The G—T 
adjustment considerably refines the probability distribution of 
the G parameter, reorienting the original data skewed toward 
lower values (owing to most of the observation data being 
collected in the colder seasons, causing an overall dip below 
the actual annual average) more toward the central sector. 

The gain curves for the remaining 18 beams are presented in 
Figure 12 (with G—T adjustment), the blue shadowed region is 
lø dispersion of G. For the full-beam calibration observations, 
a total of 17 sessions were conducted, while observations for 
M02, M05, MO8, and M14 were carried out 84 times in total 
(utilizing the DriftWithAngle mode as described in Section 2.2, 
allowing calibration of five beams in a single observation). 
Consequently, the data dispersion for these four beams is 
relatively low. Bands heavily contaminated by strong RFI are 
shaded in gray. 
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Figure 12. Corrected Gain parameter by G-T relation for outer circle beams (Beam 2 to Beam 19). The blue shadowed region indicates 1ø range. The gray shadows 


correspond to the same region as in the previous figure. 
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Figure 13. T.a difference of Beam 1 at a frequency higher than 1300 MHz, 
different lines show the percentage of the difference between each T.a file and 
Tea, on 2019 January 15. Tocai is smoothed with the frequency range of 50 MHz. 
The difference of G is shown in a red line with square markers. 


5. Discussions 


Given the measured G-T relationship, the differences 
between G and G29 become intricate. It is further 
complicated by the fact that the measurements of Gy929 for 
various beams were conducted on different dates. Mismatches 
exist in both the shape of the gain curve and the values of G. 
Within this section, we first discuss the possible origins of the 
G-T relationship, then discuss other origins of absolute flux 
calibration, and then present our online database for the Gain 
calibration. More nuanced discussions regarding result com- 
parisons are presented in Appendix. 


5.1. Possible Origin of the G-T Relation 


Various factors can contribute to the fluctuations of G, such 
as changes in the telescope response, instability of the noise 
diodes, and pointing errors. However, as shown in Section 2.3, 
the pointing is unlikely to cause significant errors in calibrator 
observation. Hence, the origin of the G—T relationship obtained 
in Section 4.2 becomes intriguing. 


5.2. Hardware Update During the Observing Epochs 


Variations in the shape of the gain-curve profiles even go 
beyond the adjustments made through the G—T relationship 
correction above in Figure 10. We find that these discrepancies 
could be a result of inaccuracies introduced by the noise diode 
files during the calibration process of Tea as discussed in 
Section 3. Differences in the shape of the gain curve may be 
due to the different conversion coefficients of the noise diodes 
(Teca) used in Equation (3). The variations in Tsa among three 
different versions (20190115, 20200531, and 20201014) are 
illustrated in Figure 13. The first version is used to determine 
Gnozo and the others are used in our data processing. We 
observed that the differences in Teaı (shown by the blue lines) 


10 


Liu et al. 


closely correspond to the differences in G indicated by the red 
line in Figure 13, in most frequency ranges, except at 1150 
MHz, where serious interference from RFI is presented. Taking 
into account the numerous sources of error in the flux 
calibration process, we suggest that the disparities between 
Gj2020 and our calculated gain curve primarily arise from the 
adjustments made through the G—T relationship correction and 
the calibration of the noise diode, although also being 
influenced by the presence of irregular RFI signals. 


5.3. Temperature-sensitive Telescope Responses 


However, the influence exerted by the Tea files tends to remain 
consistent over a given period and does not typically manifest as 
periodic effects. Notably, the relative deviations observed in the 
more recent two versions of T.a files depicted in Figure 10 
consistently exhibit the same directional bias. Therefore, it is 
improbable that they would result in scenarios where one aspect 
of G increases while another diminishes. Consequently, we 
hypothesize that the primary source of the G-T relationship may 
derive more prominently from the inherent response character- 
istics of the telescope itself. It is challenging to pinpoint precisely 
which part of the telescope response is affected by this influence. 
Given the stringent correlation between G and atmospheric 
temperature, we speculate that temperature-induced thermal 
expansion and contraction of certain hardware structures within 
the telescope might introduce interference during observations, 
such as defocus or reduced panel reflection efficiency, ultimately 
manifesting as periodic fluctuations in G. Detailed hardware 
examination is beyond the scope of this work. 


5.4. Intrinsic Flux Variation from Flux Calibrators 


Additionally, according to the flux density calibration 
guidance for the VLA,’ the calibrator 3C48 has been under- 
going a flare since 2018 or so. Other instruments beyond VLA 
have indicated measurements of the flare’s impact intensity at 
certain frequencies, with the L-band showing a relatively minor 
5% influence (NARO 2021). The flux variations induced by the 
flares of 3C48 may also be reflected in long-term observations 
of G. The blue dashed line in Figure 9 represents the 
discretization of the G values over extended observations after 
eliminating the influence of the G-T relationship, with this 
value tending toward ~5%. However, this discretization should 
already encompass both the fluctuations in the flux density of 
3C48 itself and the systematic gain errors during FAST 
observations. Thus, we infer that the impact of the flares of 
3C48 on flux density calibration is expected to be less than 5% 
and unlikely to significantly affect our other findings. 


> https: //science.nrao.edu /facilities /vla/docs /manuals /oss /performance / 
fdscale 
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5.5. Guideline of Utilizing the Detailed Data 


Detailed data on the G values for 19 beams covering the 
frequency range 1000-1500 MHz can be found on the 
HiFAST homepage: https: //hifast.readthedocs.io/fluxgain. 

These data are suitable for direct application in flux density 
calibration of observation targets with a zenith angle near 
10.15° and an atmospheric temperature close to 15° Celsius. 
When the zenith angle differs from this orientation, conversion 
using the data and formulas provided in Jiang et al. (2020) is 
necessary. Furthermore, significant deviations from the stan- 
dard temperature require adjustment for the G-T relationship 
using Equation (7). Access to atmospheric temperature data 
during observations can be obtained by contacting the FAST 
team. Note that Equation (6), fitted to obtain G, is applicable 
only for observations requiring flux density calibration 
conducted during the time depicted in Figure 6. 


6. Conclusions 


In this study, we conducted a comprehensive investigation 
into the performance of the flux density gain parameters (G) of 
FAST in different observation modes and checked the variation 
of gain in long-term observations. Our findings provide 
valuable insights into the error range of parameter G in data 
processing. We find that the differences in G obtained by 
different observation modes are negligible. Also, the variation 
of the gain parameter is small even over a timescale of years, 
although the presence of strong RFI will reduce the reliability 
of G significantly, e.g., in the frequency range from 1150 to 
1300 MHz. 

Furthermore, our analysis of the long-term fluctuation of G 
revealed a linear correlation between G and atmospheric 
temperature. We showed that correcting for temperature- 
induced errors in G can substantially reduce the dispersion of 
G in long-term observations. After this correction, the lo 
variance of G could be reduced to less than 5% for all beams. 
The gain curves for all 19 beams within the frequency range of 
1000-1500 MHz have been provided and are available for 
download from the Hi FAST homepage. 

Additionally, the test presented in Appendix B behooves us 
to note that when engaging in data processing, it is imperative 
to utilize the noise diode files temporally closest to the 
observation time. Any deviation from this established best 
practice could potentially induce errors, thus compromising the 
reliability and validity of our findings. 
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Appendix 
Comparison with G 72020 


In this section, we delve into a detailed comparative analysis 
with G2020, focusing on the implications of correcting the G-T 
relationship on beams surrounding the central beam. Figure A2 
illustrates the gain curve of the uncorrected peripheral beams, 
with the 1g and 30 ranges shown as shaded blue regions. The 
red data points, along with the error bars, represent the 
calculated Gjo929 values at the corresponding zenith angles. 
Our findings are mainly consistent with Gj2929 within one o 
range for most beams. However, beams 2, 7, 15, and 16 
indicate values exceeding Gj2929 by more than one sigma in 
certain frequency ranges. It is important to note that the 
parameters G for different beams based on Gj2929 are obtained 
from observations on various dates, making them susceptible to 
fluctuations in environmental temperature, thereby resulting in 
significant overall deviations. Figure A1 presents the rectified 
data corresponding to the parameters G for beams 2, 7, 15, and 
16 after adjustment of the G—T relationship. Subsequently, to 
this adjustment, our observational results demonstrate a closer 
alignment with those of Gj2920, with most observational points 
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falling within the lo error margin. It should be noted that we 
cannot get the temperature of the exact time for the observation 
carried out in Gy2920, the temperature data from observations 
taken around the same time in the following year are derived 
from non-systematic environmental temperature recordings at 
that specific point in time. 
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